We analyzed the C i lines associated with the damped 
Lya system observed at Zabs = 1-15 in the spectrum of 
HE 0515-4414 to derive the ^"^G/^^C ratio. The spec- 
trum was obtained by means of the UV- Visual Echelle 
Spectrograph (UVES) at the ESO Very Large Telescope 
(VLT). The obtained lower hmit ^"^G/^^C > 80 (2cr C.L.) 
shows for the first time that the abundance of ^^C in the 
extragalactic intervening clouds is very low. This rules out 
a significant contribution from intermediate-mass stars to 
the chemical evolution of matter sampled by this line of 
sight. The estimated low amount of ^^C is in agreement 
with low abundances of nitrogen observed in damped Lya 
systems - the element produced in the same nuclear cycles 
\^ • and from about the same stars as ^'^C. 
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Abstract. 



1. Introduction 

The evolution of the chemical composition of matter in the 
Universe is closely related to the history of star forma- 
tion and destruction, stellar nucleosynthesis and enrich- 
ment of the interstellar/intergalactic medium (ISM/IGM) 
with processed material. The behavior of the stable CNO 
isotopes during the cosmic time is of particular interest 
since their ratios trace the production of primary and sec- 
ondary elements^ which depends on the principal chemi- 
cal evolution parameters: the stellar initial mass function 
(IMF), the rate of the mass loss from evolved stars, and 
the dredge-up episodes leading to the mixture of the sur- 
face layers with deeper layers of the star (Wannier 1980; 
Renzini fc Voli 1981; Marigo 2001; Meynet et al. 2005). 

"^^C is a primary product of stellar nucleosynthesis and 
is formed in the triple-a process during hydrostatic he- 
lium burning (van den Hoek & Groenewegen 1997; El Eid 
2005). -^^G is supposed to have mainly a secondary origin, 
and is produced in the hydrogen burning shell when the 
GN cycle converts ^^G into ^^G (Wannier 1980). However, 
chemical evolution models considered by Prantzos et al. 
(1996) match observations better if a mixture of pri- 
mary-|-secondary origin in intermediate-mass stars is as- 
sumed for ^■'G. Evolution of massive rotating stars at 
very low metallicities may also contribute to primary ^^G 
(Meynet et al. 2005). The predicted values for the ra- 
tio TZ =^2G/^^G in the rotationally enhanced winds di- 
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* Based on observations performed at the VLT Kueyen tele- 
scope (ESO, Paranal, Chile). The data are retrieved from the 
ESO/ST-ECF Science Archive Facility 

^ A chemical element is called primary /secondary if its mass 
fraction ejected from a star is insensitive/ sensitive to the orig- 
inal metallicity of the star (Talbot & Arnett 1974). 



luted with the supernova ejects are between 100 and 4000, 

whereas the '^'^G yields of massive non-rotating stars are 
negligible, 7?. ~ 3 x 10^. Besides, Meynet et al. show that 
the rotating AGB and massive stars have about the same 
effects on the isotope production. The difference between 
them is only in the isotope composition of the massive star 
wind material and the AGB star envelopes: the former is 
characterized by very low values of 7?. ^ 3 — 5, while the 
latter have TZ ranging between 19 and 2500 (the lower TZ 
values correspond to the most massive AGB stars). 

Ghemical evolution models predict a decrease of the 
isotope ratio ^^G/^^G with time and an increase with 
galactocentric distance at a fixed time (Audouze et 
al. 1975; Dearborn et al. 1978; Tosi 1982; Romano & 
Matteucci 2003). 

For instance, a photospheric solar ratio 7^ = 95 ± 5 
(Asplund et al. 2005), representative of the local ISM 5 bil- 
lion years ago, is higher than the present value 7?, ^ 60 — 70 
obtained through optical, UV, and IR absorption line ob- 
servations as well as radio emission line measurements 
(see, e.g., Hawkins & Jura 1987; Genturion & Vladilo 
1991; Genturion et al. 1995; Goto et al. 2003; and refer- 
ences therein). It should be noted, however, that on scales 
of ^ 100 pc the local ISM is probably chemically inho- 
mogeneous (Gasassus et al. 2005) making this difference 
uncertain. 

The high resolution observations {X/6X > 50 000) of 
quasar absorption-line spectra available nowadays at large 
telescopes allow us to probe the isotope composition of 
the intervening damped Lya (DLA) systems through the 
analysis of G l lines, and, hence, to perform an important 
test of models of stellar nucleosynthesis outside the Milky 
Way. The most convenient C i transitions are from the 2u 
(A1657 A) and 3u (A1561 A) multiplets (Morton 2003) 
which cover the redshift range from 2;'~lto2;~3.5 
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Fig. 1. Normalized intensities (dots with Icr error bars) 
of the Cl Unas selected from the VLT/UVES spectrum 
of HE 0515-4414. The zero radial velocity is fixed at 
z = 1.150789 (QBR). The over-plotted smooth curves 
show the best fitted synthetic profiles {Xmin = 1-250, the 
number of degrees of freedom v = 26). The dashed verti- 
cal lines mark positions of the C i sub-components. Note 
that the A1560 absorption feature consists of two sub- 
components with different oscillator strengths: /1560.68 = 
0.0581 and /1560.70 = 0.0193. The ^^Ci* oscillator 
strengths, laboratory vacuum wavelengths and isotopic 
shifts (IS) are taken from Morton (2003) 

(look-back time ~ 7.7 — 11.8 Gyr)^ in optical spectra of 
distant QSOs. 

The analysis of the isotopic composition of the DLA 
systems is also connected to the interpretation of the hy- 
pothetical variation of the fine-structure constant (a = 
/he) at early cosmological epochs (Ashcnfelter et al. 
2004a,b; Fenner et al. 2005, hereafter FMG). 

In the present letter we report constraints on the car- 
bon isotope abundance in the Zabs =1-15 sub-DLA system 
toward HE 0515-4414 (Reimers et al. 1998). The metalhc- 
ity of this sub-DLA is found to be lower than solar: [Zn/H] 
~ -0.993 according to de la Varga et al. (2000), or [Zn/H] 
~ —0.49 as measured by Reimers et al. (2003) who revised 
the H I column density. Sects. 2 and 3 describe the obser- 
vations and the analysis. The obtained results and their 
astrophysical implications are discussed in Sect. 4. 

2. Observations and data reduction 

The observations were acquired with the UV- Visual 
Echcllc Spectrograph (UVES) at the VLT 8.2 m telescope 
at Paranal, Chile, and the spectral data were retrieved 
from the ESO archive. The seven selected exposures were 
described in detail in Levshakov et al. (2005a, hereafter 
Paper I). 

The C I lines analyzed in the present paper are located 

near the center of the echcUe orders. This minimizes possi- 
ble distortions of the line profiles caused by the decreasing 
spectral sensitivity at the edges of echelle orders. 

^ Assuming the cosmological expansion parameterized by 
to = 13.8 Gyr, h = 0.7, fi„ = 0.3, and = 0.7. 

3 [A/B] =log(iVA/iVB)-log(iVA/iVB)0. 



Following Paper I, where all details of the data reduc- 
tion are given, we work with the reduced spectra with the 
original pixel sizes in wavelength. The residuals of the cali- 
brations give (Trms ^ 1 T^^- The observed wavelength scale 
of each spectrum was transformed into vacuum, heliocen- 
tric wavelength scale (Edlcn 1966). The instrumental pro- 
file is dominated in this case by the slit width (0.8 arcsec). 
The spectral resolution calculated from the narrow lines 
of the arc spectra is FWHM = 5.60 ± 0.10 km s"! . 

After the normalization to the local continuum, the 
spectra from different exposures were rebinned with the 
step equal to the mean pixel size. To eliminate instrumen- 
tal velocity calibration errors, we used the first exposure 
as a reference frame and calculated residual differences in 
the radial velocity offsets of other exposures through the 
cross-correlation analysis. The aligned C i spectra were co- 
added with weights inversely proportional to (Tcota- The 
resulting average C I profiles are shown in Fig. 1 by dots 
with 1(7 error bars. 



3. Analysis 

The key points in the study of the isotopic abundance 
from the UV C l lines are as follows: 

1. At high redshift, Ci is observed in the most dense and 
cool sub-components of DLA systems where H2 lines are 
usually also detected. The kinetic temperature in such C i- 
bearing clouds is low. For our particular case Tkin ~ 200 K 
(Reimers et al. 2003). It means that the thermal width of 
the carbon lines is small, 6th ~ 0.5 km s~^ , and the line 
profiles can be very narrow if the turbulent broadening is 
not significant. 

2. The isotopic shift (IS)'' between and lines 
shows different values and signs (Morton 2003). For 
instance, IS1657.907/1657.916 — 1-61 km s~^ , whereas 
ISi560.682/i560.664 = -3.46 km s'^ . Although with the 
spectral resolution of FWHM = 5.6 km s~^ these lines can- 
not be completely resolved, the presence of ^"'C can shift 
and broaden ^^C lines leading to the disparity of line po- 
sitions from different multiplets. The effect is small when 
C I lines arc optically thin. The centroid of a pair of lines 
i2C-|-"C is shifted by Av = IS/(l-f 7^) with respect to the 
position of the ^^C line. Suggesting IZ = 60, this gives for 
the 1657.907 and 1560.682 lines Au = +0.026 and -0.057 
km s~^ , making the disparity of line positions compared 
to the assumption of pure ^^C equal to 0.083 km s~^ . 
However, the disparity increases with increasing optical 
depth. For saturated lines with to{}^C) > 1 and a pure 
rectangular velocity profile, the short wavelength side of 
this box-shaped absorption can be displaced by +IS or 
-IS. The centroid of ^^C+'^^G is then shifted by half of 
IS, irrespective of IZ. In this case the disparity is equal 
to 2.535 km s~^ , which already detectable at our resolu- 
tion. Note that for HE 0515-4414 the line centering can be 
done with an accuracy better than 1/10 pixel size which 



Defined as IS = c(Ai3 — Ai2)/Ai2 (in km s ^ ). 
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Fig. 2. Confidence intervals in the plane calcu- 

lated from the simultaneous fit of the C l* lines shown in 
Fig. 1. 



is about 150 m s~^ (Paper I). Since the strengths of car- 
bon lines from the UV multiplcts cover a wide dynamical 
range, the saturation effects can be correctly taken into 
account. 

3. Since the atomic carbon isotopic abundance seems 
to be unaffected by depletion, the measured ratio 
j-eflects the carbon isotopic ratio at a given red- 
shift. 

The VLT/UVES spectnmi of HE 0515-4414 reveals 
multi-component Cl profiles at ^abs = 1-15 (Fig. 1). Their 
diagnostic was performed by de la Varga et al. (2000) and 
later by Quast et al. (2002, hereafter QBR). All Ci lines 
are well described by a two-component Voigt profile model 
(see Fig. 1 in QBR). Besides, QBR found that the sub- 
components of the C I lines arising from the J = 0, 1, and 
2 levels of the ground state have different relative strengths 
indicating inhomogeneous physical conditions in the Cl- 
bearing cloud. 

The heterogeneity of the absorbing medium may cause 
small relative velocity shifts between C l transitions with 
different J similar to the shifts observed between H2 lines 
in the direction of C Ori A (Jenkins & Peimbert 1977). 
This would add an additional noise in the measurements 
of the disparity of the Ci line positions. To avoid this 
putative effect we selected from the observed eight Ci 
lines (QBR) only two J = 1 transitions which show most 
pronounced isotopic shifts. These are just the lines C i* 
A1560.68 and A1657.90 mentioned above. Other C i* lines 
(A1656.26 and A1657.37), as well as the lines with J = 
and 2 are less sensitive to the presence of "'^^C since their 
isotopic shifts are smaller than those of A1560.68 and 
A1657.90. 

Following QBR, we use two-component model to con- 
strain the isotope ratio from the minimization of the 
objective function defined by eq.(2) in Paper I. The x^j^ 
values (normalized per degree of freedom) are calculated 
in the interval 30 < 7^ < 3000 and shown as a function 
of 72. in Fig. 2. This function gradually decreases with in- 



creasing Tl and tends to a limit Xum = 1-250 at 7?. S> 1. 
This global minimum xfim = 1-250 lies within \a uncer- 
tainty range since at = 26 the expected mean value of 
X^ is equal to 1 ± 0.277. The optimized values and for- 
mal uncertainties of the model parameters for both Ci* 
absorption components (labeled by subscripts '1' and '2') 
at zero ^'^C abundance are as follows: the column density 
TVi = (3.5 ± 0.3) X 10^3 cm-2 , TVa = (6.3 ± 0.4) x lO^^ 
cm~^ , the Doppler parameter h\ = 1.10 ± 0.06 km s~^ , 
62 = 3.40 ±0.17 km s^^ , and the radial velocity difference 
At;2-i = 9.23 ±0.28 km s~^ . The corresponding synthetic 
profiles are shown in Fig. 1 by smooth curves. The model 
parameters do not change for TZ > 300. 

The dashed horizontal lines in Fig. 2 mark the la and 
2cr confidence levels. Since the curve Xmin(''^) almost flat 
at 72 > 300, the present data allow us to constrain only 
the lower limits on the carbon isotopic ratio: TZ > 235 
(1(7 C.L.), or a more conservative restriction 72 > 80 (2(t 
C.L.). The former value is, however, less certain due to a 
small gradient of the Xmin(^) function at large 72. 

We also used a 3-component model to test robust- 
ness of the obtained results. Being applied to the lines 
of Ci AA1560.3, 1658.9, Ci* AA1560.7, 1657.9, and Ci** 
AA1561.3, 1657.0, 1658.1 this model yields 72i^ > 145 and 
722a > 78 with Xmin = 1-020 {ly = 93) at 72 > 1. The 
concordance of all available C i profiles is demonstrated in 
Fig. 3. 



4. Discussion and Conclusions 

We re-analyzed the profiles of the C i lines associated with 
the sub-DLA system observed at Zabs = 1-15 in the spec- 
trum of HE 0515-4414. Our main purpose was to set a 
constraint to the carbon isotopic abundance outside the 
Milky Way, in distant intervening clouds at the cosmologi- 
cal epoch corresponding to 8.2 Gyr of the look-back time. 
Similar tasks were discussed by Carlsson et al. (1995), 
Labazan et al. (2005), and FMG. 

The present Ax^ analysis gives ^^C/^^C > 80 {2a 
C.L.). This low abundance of ^"^C does not support the 
enrichment of gas by the wind from rotating massive stars 
(Meynet et al. 2005). Besides, it can constraint the frac- 
tion of intermediate-mass (IM) stars (4 < M/Mq < 8) in 
the IMF. For instance, the FMG models for the chemical 
evolution of the Galaxy with normal IMF attain the solar 
ratio i2c/i-^C ~ 90 already at [Fe/H] « -2 for the outer 
radius and [Fe/H] « —1.5 for the solar radius case, re- 
spectively. On the contrary, the models with an enhanced 
population of IM stars produce oversolar abundances of 
13c at all metalhcities [Fe/H] < -0.3 (cf. Fig. 9 in FMG). 
Thus, these models can be ruled out by our measurements. 

Closely related to ^^C/^'^C ratio is the problem of ni- 
trogen, i^N, which is expected to be produced in the same 
way and from about the same stars as ^^C. Measurements 
in DLAs show very low N abundances which is consis- 
tent with the absence of an enhanced population of IM 
stars to the chemical evolution of the gas (Centurion et al. 
2003). Thus, low ^^C and low ^^N seem to be in agreement. 
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It should be noted, however, that the chemical evolution 
of nitrogen is still not well understood since the existing 
models predict abundances higher than observed in DLAs 
(see, e.g.. Fig. 6 in FMG). 

The present bound to the abundance of ^^C and, as a 
result, to the contribution of the AGB stars to the chem- 
ical evolution of the Zabs — 1-15 system may be indirectly 
related to the claims of a possible time variation of the 
fine-structure constant, a (see Murphy et al. 2004 and ref- 
erences therein). A significant portion of the sample used 
by Murphy et al. involves the comparison of Mg 11 and Fe 11 
wavelength shifts. Later on, it was shown that over-solar 
abundances (^ 0.3 dex) of ^^'^^Mg isotopes with respect 
to ^^Mg in the absorbing material can imitate an appar- 
ent variation of a in the redshift range between 0.5 and 
1.8 (Ashenfelter et al. 2004a,b). The production of Mg iso- 
topes is believed to occur in about the same stars which 
produce ^^C. Thus, the obtained bound to the amount of 
^■^C, taken as a typical value for the DLA systems^, poses 
a limit to the role of the AGB stars in mimicking the a 
variations. However, the variability of a has not been sup- 
ported by more recent studies (Quast et al. 2004; Chand et 
al. 2004; Levshakov et al. 2005a, b). This may suggest the 
presence of other systematic errors to explain the different 
results concerning the variability of a. 
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Fig. 3. Same as Fig. 1 but for a 3-componcnt Voigt profile model fitted to all available C i lines. The over-plotted 
smooth curves show the best fitted synthetic profiles (Xmin = 1-036, the number of degrees of freedom v = 125). The 
dashed vertical lines mark positions of the C i sub-components 



